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ABSTRACT 

Context. Some merging galaxy clusters host diffuse elongated radio sources, also called radio relics. It is proposed that these radio 
relics trace shock waves in the intracluster medium (ICM) created during a cluster merger event. Within the shock waves particles 
are accelerated to relativistic energies, and in the presence of a magnetic field synchrotron radiation will be emitted. Here we present 
Giant Metrewave Radio Telescope (GMRT) and Westerbork Synthesis Radio Telescope (WSRT) observations of a new double relic 
in the galaxy cluster ZwCI 0008.8+5215. 

Aims. The aim of the observation is to understand the phenomenon of radio relics. 

Methods. We carried out radio continuum observations at 241 and 610 MHz with the GMRT, and 1.3 - 1.8 GHz observations with 
the WSRT in full polarization mode. Optical V, R, and I band images of the cluster were taken with the 2.5m Isaac Newton Telescope 
(INT). An optical spectrum, to determine the redshift of the cluster, was taken with the William Herschel Telescope (WHT). 
Results. Our observations show the presence of a double radio relic in the galaxy cluster ZwCI 0008.8+5215, for which we find a 
spectroscopic redshift of z = 0.1032 + 0.0018 from an optical spectrum of one of the cD galaxies. The spectral index of the two relics 
steepens inwards to the cluster center. For part of the relics, we measure a polarization fraction in the range ~ 5 - 25%. A ROSAT 
X-ray image displays an elongated ICM and the large-scale distribution of galaxies reveals two cluster cores, all pointing towards 
a binary cluster merger event. The radio relics are located symmetrically with respect to the X-ray center of the cluster, along the 
proposed merger axis. The relics have a linear extent of 1.4 Mpc and 290 kpc. This factor of five difference in linear size is unlike that 
of previously known double relic systems, for which the sizes do not differ by more than a factor of two. 

Conclusions. We conclude that the double relics in ZwCI 0008.8+5215 are best explained by two outward moving shock waves in 
which particles are (re) accelerated trough the diffusive shock acceleration (DSA) mechanism. 

Key words. Radio Continuum: galaxies - Galaxies: active - Clusters: individual : ZwCI 0008.8+5215 - Cosmology: large-scale 
structure of Universe 

1, Introduction Some smaller radio relics (< 500kpc) have been been ex- 

plained by old ("fossil") radio plasma from a previous episode 
of AGN activity. T hese sources are called AGN relics (see 
> - Radio relics are filamentary structures often located in the iKempner et alJl2004j, for an overview of the classification of dif- 
•Sh periphery of galaxy clusters. It is proposed that large ra- fuse radio sources). The fossil radio pla sma could also have been 
^ dio rel ics trace shock w a ves generated by cluster merger compressed, creating a radio phoenix (EnBlin & Gopal-Krishna| 
Jh ■ events (lEnsslin et all fl998l: iMiniati et al.l [200lh . At the shock l2QQll l EnBlin & Bruggen| | 2002|) . Both radio phoenices and AGN 
front particles from the thermal gas are accelerated to rel- rellcs > are characterized by a very steep (a < -1.5, F v oc v a , 
ativistic energies by DSA mechanism i n a fi rst-ord er Ferm i where a 1S the spectral index) and curved radio spectra due to 
process f Krymskii \l91% [A xford et^D fl977l: 1^ [197&K synchrotron and Inverse Compton (IC) losses 
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Blandford & Ostrikerl | 1978[ iDrurvl Il983i iBla ndford & Eic hlerl In the hierarchical model of structure formation galaxy clus- 
19871: U ones & Ellison 1991; Malk ov & O'C Drurvll2001h T In the ter grow by the accretion of gas from the surrounding interga lac- 



presence of a magnetic field these particles emit synchrotron ra- tic medium (IGM) and through mergers with other clusters and 

diation at radio wavelen gths. Another possibility mentioned by galaxy groups. Large radio relics are exclusively found in dis- 

iMarkevitch et al.l (|2005), is that the shock re-accelerates rela- turbed clusters, indicative of merger activity. This supports the 

tivistic fossil electrons injected previously into the ICM by for idea that shocks generated during cluster merger events can be 

example AGN. They note that from an observational point of responsible for the non-thermal radio emission. Hydrodynamical 

view, this case will probably be indistinguishable from the direct models of structur e formation, including particle acceleration 



shock acceleration mentioned above. An alternative scenario has mechanisms (e.g., iHoeft & Briiggenl 120071: iHoeft et all 2008; 
been recently proposed, namely th at relics arise from emission iPfrommen 1 120081: iBattaglia et al.l 12009) make predictions about 
of secondary cosmic ray electrons (Keshet 20 10|). the location, orientation and radio power of relics in merging 
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cluster s. Amongst the several dozen radio relics known to date 
(e.g., [va n Weere net al.ll2009bc iRudnick & L emmerman 2009; 
Pizzo & de Bruvn 120091: Giacintu cci et al.l 120081: [B rentiens 
2008; Kempner & SarazidBOOl. iGiovannini et al.lll999l Il99ll) . 
there are a few rare double relic systems, with two relics lo- 
cated symmetrically on opposite sites of the cluster center 
(e.g., van Weeren et al. 2010; Brown et al. 2010; Bonafede et al. 
2009; va n Weeren et al ll2009aUVenturi et al.ll2QQ7b iBagchi et al.l 
20061 iRottgering et al.lfl997h . These double relics can be used 
to constrain the mer ger geometry and timescales involved 
(iRoettiger et al.l fl999l) . If radio relics trace outward traveling 
shock waves in which DSA takes place, then the radio plasma 
in the post-shock region should have a steeper spectrum due 
to IC and synchrotron losses. For a relic seen close to edge- 
on, the luminosity profile across the width of the relic can 
then be used to constr ain the magnetic fields strength at the lo- 
cation of the shock (iMarkevitch et al.l 120051: iFinoguenov et al.l 
l2oTob , This is because the downstream luminosity profile should 
directly reflect the synchrotron l osses, which in turn dep end 
on the magnetic field strength. In Ivan Weeren et al.l (120101) we 
presented observations of a double relic in the merging clus- 
ter CIZA J2242.8+5301, which provided evidence for DSA in 
galaxy cluster merger shocks. For the largest relic we derived 
a magnetic field strength of about 5-7 //Gauss by modeling the 
relic's luminosity profile across the width of the relic (although a 
strength of about 1.2 //Gauss could not be completely ruled out). 

Because there are only a few double relics systems 
known, we carried out an extensive search in the 1.4 GHz 
NVSS (ICondon et al.ll998l) . 32 5 MHz WENSS dReneelink et al.1 
119971) . and 74 MHz VLSS jCohen et al.1 12007|) surveys for 
the presence of arc-like radio structures around X-ray selected 
galaxy clusters from the ROSAT All-Sky Survey. This search al- 
ready resulted in t he disc overy of several new radio relics, see 
Ivan Weeren et al . (20093513).' In this paper we present the dis- 
covery of a double relic in the galaxy cluster ZwCl 0008.8+5215 
which showed faint elongated structures in both NVSS and 
WENSS images. 

The layout of this paper is as follows. In Sect. [2] we give an 
overview of the observations and the data reduction. In Sect. [3] 
we present the radio and spectral index maps as well as opti- 
cal images around the radio sources. In Sect. [4] we discuss the 
merger scenario and the magnetic field strength at the location 
of the relics. We end with conclusions in Sect. [5] 

Throughout this paper we assume a ACDM cosmology with 
Hq = 71 km s" 1 Mpc -1 , £_ m = 0.3, and Q\ = 0.7. All images are 
in the J2000 coordinate system. 

2. Observations & Data Reduction 

2.1. GMRT observations 

The GMRT observations were taken in dual-frequency mode, 
recording RR polarization at 610 MHz and LL polarization at 
241 MHz. Total useable bandwidth was 30 MHz at 610 MHz 
and 6 MHz at 241 MHz. The GMRT software backend (GSB; 
iRov et al.1 2009) was used giving 512 spectral channels. The to- 
tal on source time was 220 min (3.7 hr). A summary of the ob- 
servations is given in Table [T] 

For the reduction we used the NRAO Astronomical Image 
Processing System (AIPS) package. The 610 MHz dataset was 
visually inspected for the presence of radio frequency interfer- 
ence (RFI) which was subsequently removed. Strong RFI on 
short baselines was present in the 241 MHz band. This RFI was 
fitted and subtracted from the data using the technique described 



Table 1. GMRT observations 





241 MHz 


610 MHz 


Observation date 


JNovember zz, zUU^ 


JNovember zz, zUUy 


Usable bandwidth 


6 MHz 


30 MHz 


Channel width 


62.5 kHz 


62.5 kHz 


Polarization 


LL 


RR 


Integration time 


8 sec 


8 sec 


Total on- source time 


220 min 


220 min 


Beam size 


14.65" x 12.5" 


6.3" x 5.3" 


Rms noise (cr rms ) 


483 / Jy beam -1 


38 yL/Jy beam -1 



byjAthreval ([2009) which was implemented in Obit (iCottonl 
120081) . Summarized, the fringe- stopped correlator output of a 
baseline oscillates with the fringe-stop period in the presence of 
RFI. The fittin g routine tries to fit such a signal and subtracts it 
from the data (lAthreyall2009l) . This has the advantage that most 
of the visibilities from the short baselines are not removed (i.e., 
flagged) and hence spatial structure on large angular scales is 
preserved. We assumed the RFI to be constant over periods of 
10 min and required a minimal amplitude of 5 Jy for the RFI 
to be subtracted. Further visual inspection of the data was car- 
ried out to remove some remaining RFI which could not be fitted 
because of short timescale variations of the RFI signal. 

Amplitude and phase corrections were obtained for the cali- 
brator sources using 5 neighboring channels free of RFI. These 
solutions were applied to the data before determining the band- 
pass. The bandpass was then applied and gain solutions for the 
full channel range were determined . The fluxes of the prim ary 
calibrators were set according to the lPerley & Taylor! (Il999l) ex- 
tension to the lBaars et al.l (1 1977b scale. Several rounds of phase 
self-calibration and two final rounds of amplitude and phase self- 
calibration were carried out. We used the polyhedron method 
dPerleyll 9891: ICorn well & Per ley 1992) for making the images to 
minimize the effects of non-coplanar baselines. Both in the 610 
and 241 MHz images, there were several bright sources in the 
field that limited the dynamic range. Direction dependent gain 
solutions for these sources were obtained and these sources were 
subtracted from t he data. This me thod is commonly referred to 
as "peeling" (e.g jNoor dam 2004]). 

Final image s were made using "briggs" weighting (with ro- 
bust set to 0.5. lBriggslll995l) . Images were cleaned down to 2 
times the rms noise level (20-^) within the clean boxes. The 
images were corrected for the primary beam response^. The un- 
certainty in t he calibration of the a bsolute flux-scale is between 
5 - 10%, see lChandra etaD (12004 . 

2.2. WSRT 1.3-1.7 GHz observations 

WSRT observations of ZwCl 0008.8+5215 were taken in the L- 
band, see TableO The observations were spread out over various 
runs of several hours, resulting in more or less full 12 hours syn- 
thesis coverage. Frequency switching, between the 18 and 21 cm 
setups, was employed every 5 min to increase the spectral base- 
line of the observations. Each setup had a total bandwidth of 
160 MHz, evenly divided over 8 sidebands (IF). All four polar- 
ization products were recored with 64 channels per IF. 

The data were calibrated using the CASA0 package. The 
L-band receivers of the WSRT antennas have linearly polar- 

1 http://grnrt.ncra.tifr.res.in/gmrt_hpage/Users/doc/manual/ 
UsersManuaVnode27.html 

2 http://casa.nrao.edu/ 
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Frequency bands 21 cm (IFs) 
Frequency bands 18 cm (IFs) 
Bandwidth per IF 
Number of channels per IF 
Channel width 
Polarization 
Observation dates 
Integration time 
Total on- source time 
Beam size 
Rms noise (cr rms ) 



1311, 1330, 1350, 1370, 1392, 1410, 1432, 1450 MHz 
1650, 1668, 1686, 1704, 1722, 1740, 1758, 1776 MHz 
20 MHz 
64 

312.5 kHz 

XX, YY, XY, XY 

28 March, 25 & 31 May, 7 & 16 June 7, 2009 
30 sec 

12.5 hr (21cm), 12.5 hr (18cm) 

23.5" x 17.0" (21cm), 23.5" x 17.0" (18cm) 

27 yt/Jy beam -1 (21cm), 33 yt/Jy beam -1 (18cm) 



ized feedfl Bandpass and gain solutions were determined from 
observations of the standard calibrators 3C48, 3C138, 3C147, 
3C286, and CTD93. Time ranges for antennas affected by shad- 
owing were removed. Some IFs were also affected by RFI and 
had to be partly flag ged. The fluxes for the calibrators were 
set accordi n g to t he iPerley & Taylor! (|T999) extension to the 
Baa rs et aD (119771) scale. Frequency (channel) dependent leak- 
age terms (D-terms) were calculated from observations of a 
bright unpolarized point source, while the polarization angles 
were set using either 3C138 or 3C286. We assumed -66.0° and 
15.0° for the R-L phase difference of 3C286 and 3C138, respec- 
tively. The data were then exported into AIPS and several rounds 
of phase self-calibration, followed by two rounds of amplitude 
and phase self-calibration, were carried out. The solutions for 
the amplitude and phase self-calibration were determined com- 
bining both XX and YY polarizations because Stokes Q is not 
necessarily zero. Separate images for each IF were made (16 in 
total) using robust weighting of 0.75 (21 cm) and 1.25 (18 cm) 
giving roughly the same resolution. These images were corrected 
for the primary beam attenuation (A(r)) 



A(r) = cos 6 (r x y x C) , 



(1) 



with r the distance from the pointing center in degrees, v the 
observing frequency in GHz and C = 68, a constanfl. 

Images were cleaned to about and clean boxes were 

used. For each frequency setup the eight images from each IF 
were combined, convolving the images to the same resolution 
of 23.5" x 17.0" and weighting the images inversely propor- 
tional to cr^ ms . The images are centered at frequencies of 1382 
and 1714 MHz. The final noise levels are 27 and 33 yuJy beam -1 
for the 21 cm and 18 cm images, respectively. 

2.3. Optical images & spectroscopy 

Optical images of the cluster were made with the Wide Field 
Camera (WFC) on the INT. The observation were taken be- 
tween 1 and 8 October, 2009. The seeing varied between 0.9 
and 1.3" during the observations. Most nights were photomet- 
ric. Total integration time was about 6000 sec per filter. The data 
were reduced in a standar d way using I RAF (iTodyll 19861 Il993l) 
and the mscred package (I Valdesl [l 99 8h . The R and I band im- 
ages were fringe corrected. The individual exposures were av- 
eraged, pixels were rejected above 30-^ to remove cosmic rays 
and other artifacts. Observations of standard stars on photomet- 
ric nights were used to calibrate the flux scale. The images have 



o'< 4x10" 



3 2X10" 1 ' 




5000 o 
wavelength [A] 



3 The WSRT records XX = I - Q, YY = I - Q, XY = -U + iV, and 
XY = -U - iV. I, Q, U, and V are the standard Stokes parameters 



4 from the WSRT Guide to Observations 



Fig.l. Spectrum of the galaxy 2MASX J00112171+5231437 
taken with the WHT ACAM V400 grating. Various absorption 
features are indicated. 



a depth (signal to noise radio (SNR) of 5 for point sources) of 
approximately 24.5, 23.9, 23.3 magnitude (Vega) in the V, R and 
I band, respectively. 

To determine the redshift of the cluster a 600 sec WHT 
ACAM long- slit spectrum, with the V400 grating, w as taken 
of 2MASX J001 12171+5231437 dSkrutskie et alJl2003b . On the 
INT images this galaxy was identified being the largest cD 
galaxy of the cluster. The spectrum was taken on November 2, 
2010 with a slit width of 1.5". Standard long-slit calibration was 
done in IRAF. 



3. Results 

3. 1 . Redshift of ZwCl 0008. 8+52 1 5 

No redshift is available in the literature for the galaxy clus- 
ter ZwCl 0008.8+5215, with coordinates RA 00 h ll m 25.6 s , 
DEC +52°3r41". The cluster is located relatively close to the 
galactic plane at a galactic latitude b = -9.86°. Galactic ex- 
tinction is 0. 812 mag in the R-ba nd and 0.111 in the K-band 
according to Schlegel et al. ( 1998). The spectrum of cD galaxy 
2MASX J001 12171+5231437 is shown in Fig. ffl From this 
spectrum we determine a redshift of 0.1032 + 0.018 for the 
galaxy, which we adopt as the redshift for ZwCl 0008.8+5215. 
2MASX J001 12171+5231437 is associated with a complex dis- 
turbed radio source, see Sect. [ 
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3.2. Thermal ICM and galaxy distribution 

ZwCl 0008.8+5215 is seen in the ROSAT All-Sky Survey 
as an east- west elongated source, see Fig. |2l and listed as 
1RXS J001 145.3+523147 (IVoges et al.l fl999h . Using the red- 
shift and the ROSAT count rate we find an X-ray luminosity 
(L x , 0.1-2.4 kev) of ~ 5 x 10 43 erg s" 1 . With the Lx , 0.1-2.4 kev- 
temperature scaling relation from Pra tt et all (12009), Table B.2 
BCES Orthogonal Fitting Method, we find a corresponding tem- 
perature of ~ 3 to 4 keV. 

We computed galaxy iso-densities from the INT im- 
ages. We first created a catalog of objects using Sextractor 
(iBertin & Arnoutsl IT996). We then removed all point-like ob- 
jects (i.e., stars) from the catalogs. To exclude galaxies not be- 
longing to the cluster we selected only galaxies with R-I and 
V-R colors within 0.15 magnitude from the average color of 
the massive elliptical cD galaxy 2MASX J001 12171+5231437 
(see Fig. [14] and Sect. 13.11) . The range of 0.15 in the colors was 
taken to maximize the contrast of the cluster with respect to the 
fore and background galaxies in the field, but not being too re- 
strictive so that a sufficient number of candidate cluster mem- 
bers was selected. The galaxy iso-density contours are shown in 
Fig. [2] The cluster shows a pronounced bimodal structure, with 
two cores separated by about 700 kpc. The cluster extends some- 
what further in the east- west direction than the X-ray emission 
from ROSAT. The cD galaxy 2MASX J001 12171+5231437 be- 
longs to the western subcluster (i.e., it is located at the center of 
the subcluster). The eastern subcluster also hosts a separate cD 
galaxy (2MASX J00121892+5233460, see Fig. [13]). Although 
we do not have a spectroscopic redshift for this galaxy, the (i) 
color and (ii) R and K magnitudes are in agreement with a sub- 
clust er located at the same redshift as the w estern subcluster 
(e.g.. lWillottet al. 2003} |de Vries et al . 2007). The same is true 
for the other mass ive elliptical galaxies found in both subclus- 
ters, see Sect. 13.51 Therefore, both the X-ray and optical obser- 
vations point towards a bi-modal galaxy cluster, indicative of an 
ongoing merger event. As we will show in the next sections, the 
radio observations also point towards a merger scenario. 

3.3. Radio continuum maps 

The WSRT 1382 MHz image is shown in Fig. [3l It reveals a 
large arc of diffuse emission on the east side of the cluster and a 
smaller faint diffuse source on the west side of the cluster, sym- 
metrically with respect to the cluster center. We classify these 
sources as radio relics based on their location with respect to the 
cluster center, their morphology, and the lack of optical coun- 
terparts. The relics are located about 850 kpc from the center of 
the X-ray emission. Several complex tailed radio sources, related 
to AGN activity, are also visible. The WSRT 1714 MHz image 
is similar to the 1382 MHz image, although the overall signal 
to noise ratio is less, therefore revealing less of the diffuse ex- 
tended relics. The radio relics are also visible in the GMRT 610 
and 241 MHz images (Figs. [5] and©, although at 241 MHz the 
SNR on the relics is < 5 per beam. 

To facilitate the discussion we have labeled variou s sou rces 
in Figs. [3] and [7J Optical overlays can be found Sect. 13.51 The 
integrated fluxes, spectral indices, radio power and largest linear 
size for the two relics (RE & RW) are displayed in Table [3] 

Relic RE consist of two parts, a smaller region of emission 
to the north and a larger one in the south (most clearly seen 
in Figs. [5] and 0. In the 1382 MHz image the two regions are 
seen connected. The eastern boundary of RE is somewhat more 
pronounced, while on the western side the emission fades more 
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Fig. 2. X-ray emission from ROSAT, tracing the thermal ICM, is 
shown by the color image. The original image from the ROSAT 
All Sky Survey was convolved with a 200" FWHM Gaussian. 
Solid contours are from the WSRT 1382 MHz image and drawn 
at levels of [1, 2, 4, 8, . . .] x 4cr rms . The resolution of the WSRT 
image is 23.5" x 17.0". Dashed contours show the galaxy iso- 
density distribution derived from the INT images. Contours are 
drawn at [1, 1.1, 1.2, 1.3, 1.4 ...] x 0.38 galaxies arcmin -2 using 
the color cuts as described in Sect. 13.21 




Right Ascension 



Fig. 3. WSRT 1382 MHz image. Contour levels are drawn at 
[-1, 1, 2, 4, 8, . . .] x 4<x rms . Negative contours are shown by the 
dotted lines. The beam size is 23.5" x 17.0" and shown in the 
bottom left corner of the image. Sources are labeled as in Fig. 13 
The white + symbols mark the centers of the two subclusters 
based on the iso-density contours from Fig. [2] 
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Source 


S 241 MHz 


^610 MHz 


S 1382 MHz 


S 1714 MHz 


^1714 MHz 
"241MHz 


^1.4 GHz 


LLS^ 




Jy 


mJy 


mJy 


mJy 




10 24 W Hz" 1 


kpc 


RW 


0.11+0.03 


56 + 8 


11 + 1.2 


8.9 + 1.2 


-1.49 + 0.12 a 


0.37 


290 


RE 


0.82 + 0.09 


230 + 25 


56 + 3.5 


37 + 2.7 


-1.59 + 0.06 a 


1.8 


1400 


A 


0.186 + 0.020 


86 + 9.0 


37 + 1.9 


29 + 1.6 


-0.96 + 0.06 


1.1 


270 


B 


U.llz + U.UIj 


on i n c\ 
o / + 9.U 


o a i 1 n 
34 + 1.9 


OO i 1 c 
Z5 + l.O 


-U.ol + U.Uo 


a no 




C 


0.423 + 0.045 


205 + 21 


88+4.6 


69 + 3.6 


-0.94 + 0.06 


2.6 


249 


D 




2.2 + 0.3 


1.06 + 0.1 


1.2 + 0.1 


-0.63 + 0.16 c 


0.030 




E 


2.81+0.28 


1016+ 102 


387 + 20 


302 + 15 


-1.14 + 0.05 


11.9 


200 


F 


0.280 + 0.030 


135 + 14 


58 + 3 


47 + 3 


-0.93 + 0.06 


1.7 


151 


G 


0.040 + 0.005 


16 + 2 


7.3 + 0.4 


6.3 + 0.4 


-0.94 + 0.07 


0.22 


38 


H 




8.7 + 1.2 


2.6 + 0.2 


2.7 + 0.2 


-1.18 + 0.15 c 


0.081 


60 


I 




2.8 + 0.5 


1.2 + 0.19 


0.76 + 0.18 


-1.15 + 0.23 c 


0.037 


45 



a see also Fig. [15] 

b largest linear size 

c between 610 and 1714 MHz 




Right Ascension Right Ascension 



Fig. 4. WSRT 1714 MHz image. Contour levels are drawn at 
[-1, 1, 2, 4, 8, . . .] x 4<x rms . Negative contours are shown by the 
dotted lines. The beam size is 23.5" x 17.0" and shown in the 
bottom left corner of the image. 



Fig. 6. GMRT 241 MHz image. Contour levels are drawn at 
[-1, 1, 2, 4, 8, . . .] x 40-^. Negative contours are shown by the 
dotted lines. The beam size is 14.65" x 12.50" and shown in the 
bottom left corner of the image. 



slowly in the direction of the cluster center. The relic has a total 
extent of 1.4 Mpc. The surface brightness varies across the relic 
fading at the extreme northern and southern ends. The northern 
diffuse patch has a "notch" like region of higher surface bright- 
ness. Relic RW has a much smaller extent of 290 kpc. The west- 
ern boundary is more pronounced in the WSRT images. A com- 
pact source in the middle of RW is associated with a background 
galaxy. 

The cluster also hosts a number of complex radio sources 
related to AGN activity, for a short discussion on these sources 
see Sect. 13.51 



3.4. Spectral index and polarization maps 

A spectral index map was computed using both the WSRT 
and GMRT images, including only common UV ranges. Both 
the WSRT and GMRT datasets have relatively good inner UV- 
coverage. The largest detectable angular scale is limited to about 
16 r at 610 MHz, which is sufficient not to resolve out the ex- 
tended radio relics. The inclusion of maps at four different fre- 
quencies enables us to map the spectral index over the low sur- 
face brightness radio relics. Spectral index maps made with only 
two frequency images were too noisy to map the spectral index 
across the relics. The spectral index map was created by fitting 
a single power-law through the flux measurements at 241, 610, 
1382, and 1714 MHz. In this way, we only fitted for the slope 
and normalization of the radio spectrum, ensuring the number of 




Fig. 5. GMRT 610 MHz image. Contour levels are drawn at [-1,1, 2, 4, 8, ...]x 4<x rms . Negative contours are shown by the dotted 
lines. The beam size is 6.3" x 5.3" and shown in the bottom left corner of the image. 



free variables in the fit remains as low as possible (at the cost of 
detecting spectral curvature). The technique of combining maps 
at more than two frequencies has another advantage that errors 
in the maps arising from RFI, calibration errors, deconvolution 
errors, slightly different UV coverage, etc., are suppressed in the 
spectral index map as long as they do not correlate at the same 
location and spatial frequencies on the sky. 

Pixels in the spectral index map were blanked if any of cor- 
responding pixels in the individual maps fell below 1.50-^. 
Special care was taken about the precise alignment of the maps, 
we slightly shifted the GMRT maps by about a quarter of the 
synthesized beam, removing a small spectral index gradient vis- 
ible across all the point sources. The result is shown in Fig. [8] 



For relic RW, the spectral index steepens to the north and 
eastwards to the cluster center, from -0.9 to -2.0. The spectral 
index for relic RE also varies roughly between -0.9 and -2.0. 
The overall spectral index at the east side of relic RE is about 
-1.2 There is an overall trend of spectral steepening towards the 
west, see also Fig. QIJ The spectral index is correlated with the 
surface brightness of the relic, the brightest parts have a flatter 
spectral index (see also Fig.[T6l). 

The polarization map from the WSRT at 1382 MHz is shown 
in Fig. [9] No useful polarization information could be extracted 
from the WSRT 18 cm observations. The polarization map re- 
veals that most of the compact sources are polarized below the 
5% level. Some polarized emission is detected from the two ra- 
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Fig. 8. Left: Spectral index map. Spectral indices were computed by fitting a single power-law radio spectrum through the fluxes at 
241, 610, 1382, and 1714 MHz. Solid contours are from the 1382 MHz WSRT map and drawn at levels of [1, 4, 16, 64, . . .] x 6cr rms . 
The resolution of the map is 23.5" x 17.0". Right: Spectral index uncertainty map. The map is computed on the basis of <x rms values 
for the individual maps. 




Fig. 7. GMRT 610 MHz image with sources labeled. The blue 
and white + symbols mark the centers of the two subclusters 
based on the iso-density contours from Fig. 12 



dio relics, although at low SNR. For relic RW the polarization 
fraction is around 5 - 1 0% (reported polarization f ractions are 
corrected for Ricean bias: IWardle & Kronbergl 19741) . For RE the 
polarization fraction varies, with a maximum of ~ 25%. For the 
fainter parts of the relics no polarized emission is detected, but 
this is expected if the relics are polarized at the 30% level or less. 
Most polarization E- vectors are aligned perpendicular to the ma- 
jor axis of the two relics (for the parts where polarized emission 
was detected). 



Fig. 9. WSRT 1382 MHz polarization E- vector map. Total polar- 
ized intensity is shown as grayscale image. Vectors depict the po- 
larization E- vectors, their length represents the polarization frac- 
tion. The length of the E- v ectors are corrected for Ricean bias 
dWardle & K ronberg l 19741) . A reference vector for a polarization 
fraction of 100% is shown in the bottom left corner. No vectors 
were drawn for pixels with a SNR < 5 in the total polarized in- 
tensity image. Contour levels are drawn at [1 , 2, 4, 8, . . .] x So-m^ 
and are from the Stokes I 1382 MHz image. The beam size is 
23.5" x 17.0" and shown in the bottom left corner of the image. 



Spectral index and polarization properties for the compact 
sources are discussed in Sect. 13.51 
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3.5. Radio galaxies in the cluster 

The cluster hosts several interesting tailed radio sources. Radio 
overlays on optical images are shown in Figs. [TO] to [14] Sources 
are referred to as in Fig. [7] The morphology of these radio 
sources is as expected for a system undergoing a merger, with 
high galaxy velocities with respect to the ICM. The radio powers 
reported for the sou rces in Table |3] are co nsistent with them be- 
ing F R-I sources ( Fanaroff & Riley] 1 19741) located in the cluster 
fe.p. . lOwen & Whitdll99ll: lOwen & Ledlowlll994l: iBaum et al.1 

Source A is a "head-tail" source belonging to the elliptical 
galaxy 2MASX J00 120320+5234 132. The spectral index steep- 
ens along the tail, from -0.6 to -2.1 to the west. The "head" is 
polarized at the 1% level. The tail has a total extent of 270 kpc 
in the 610 MHz image. Source B is a double-lobe radio source 
(also listed as NVSS J001242+522717 or B0010.0+5210) con- 
sisting of a central core and two radio lobes. Additional diffuse 
emission is seen surrounding the lobes and the central radio core. 
The spectral index does not vary much over the core and radio 
lobes, the average spectral index is about -0.6. The radio emis- 
sion from the source has a polarization fraction between 10 and 
20%. An optical counterpart (2MASX J00124217+5227182) is 
centered on the radio core. The source could be a cluster member 
as its color is similar to that of other galaxies in the cluster. 

Source C is a narrow angle tailed (NAT) radio source 
with 2MASX J001 12859+5228096 being the optical counter- 
part. The spectral index steepens from -0.6 to -1.7 along 
the tails as expected for a radio galaxy moving eastwards 
with respect to the ICM. The radio emission is polarized 
at the 2 - 4% level. Source D belongs to the cD galaxy 
2MASX J00121892+5233460 in the eastern part of the cluster 
and has a spectral index of -0.8. No polarized emission is de- 
tected from the source. 

Source E is the brightest radio source (4C +52.01) in the 
cluster and has a very complex morphology. The radio emis- 
sion has a sharp western boundary, while the eastern boundary 
is more diffuse with "fans" of emission. The counterpart is the 
galaxy 2MASX J001 12171+5231437, which has a close com- 
panion to the west. The spectral index steepens to the north and 
south from -0.9 to -1.7. Fractional polarization for the brightest 
emission is below 1%. 

Source F has a complex morphology and is located to the 
northwest of E. There is a bridge of steep spectrum emission 

(with a 1.4) between E and F. The polarization fraction of 

F varies between 2 and 10%. The brighter parts of the source 
have a spectral index of -0.75. Because the spectral index here 
is flatter than the bridge of emission between E and F, this 
suggests the source is a separate radio galaxy with the emis- 
sion not coming from E. The most likely counterpart is the 
galaxy 2MASS J001 1 1 135+5232421. The complex morphology 
of both E and F suggest that the radio sources are significantly 
disturbed, possibly due to the merger event. The radio morphol- 
ogy of both E and F suggest the sources moves westwards with 
respect to the ICM. 

We could not identify optical counterparts for sources G, H 
and I. However, a diffraction spike in the INT images from a 
nearby bright star partly covers source G, likely blocking our 
view of the optical counterpart. The high surface brightness of G 
indicates the source is probably associated with an AGN. 



in the galaxy cluster ZwCl 0008.8+5215 




Right Ascension 

Fig. 10. Source A. GMRT 610 MHz contour are drawn at levels 
of VTH 2, 4, 8, . . Jx4cr rms . The beam size is 11.75" X7.65" and 
shown in the bottom left corner of the image. 
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Fig. 11. Source B. GMRT 610 MHz contour levels are displayed 
as in Fig.fTOl 

4. Discussion 

4.1. Origin of the double radio relic 

Hydrodynamical simulations of cluster mergers show that 
the p r ocess take s of the order of 10^ yr (e . g.. IVazza et al. 
120101: iPaul et al.1 12010b iBattaelia et all 120091: ISkillman et al. 
200a lHoeft etaLll2008l : IPfrommer et al.ll2 008: Ric ker & Sarazin 
l200ll) . During a cluster merger, "internal" shock waves are gen- 
erated. Typically, these have lower Mach numbers (M < 5) 
than "external" shock waves which are generated by the in- 
fall of unprocessed gas from the surrounding IGM and hav- 
ing M > 10. The Mach numbers of internal shocks are low 
because the sound speed in the gas of the main (bigger) clus- 
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Fig. 14. Source E & F. GMRT 610 MHz contour levels are displayed as in Fig.fTOl 



ter and the velocity of the in-falling subcluster both reflect the 
same gravitational potential of the main cluster. Merger events 
which generate shock waves with At > 3 are rare, and these 
are mainly formed in major merger events, with the mass ra- 
tio of the two merging clusters approaching unity. In the case 
of a binary cluster, merger two shocks are produced along the 
merger axis. As the shock waves propagate outward into a lower 
density environment their Mach number increases. The shock 
structure may get broken when it interacts with the galaxy fil- 
aments connected to the cluster. This may explain the pres- 
ence ofjl^^otchjike'' feature observed in the north of relic 
RE dPaul et al.ll2010l) . These notch-like str uctures are also ob- 
served for the d ouble relics in Abell 3 376 (iBagchi et al. 2006) 
and Abell 3667 (iRottgering et al.ll 19971) . 

At the location of the shock front particles are proposed 
(lEnsslin et al.| [T998) to be (re)accelerated by the DSA mecha- 
nism. In this case, the injection spectral index of the radio emis- 



sion is related to the Mach number of the shock. Behind the 
shock front particles cool through IC and synchrotron losses. 
The spectral index should thus steepen inwards to the clus- 
ter center. The overall integrated spectral is still a power-law 
though, but steeper by abou t 0.5 unit (e.g. jBlandford & Eichlerf 
119871: ISarazinl 1 19991 l2002h than the injection spectral index. 
Relics which such power-law spectra hav e been found for exam- 
ple in Abell 521 and th e Coma cluster (iGiacintucci et al.l l2008: 
iGiovannini et al.ll 19911) . 

We have fitted a power-law radio spectrum for the integrated 
fluxes of the radio relic reported in Table The fitted radio 
spectra for the two relics are displayed in Fig. [15] RE is well 
fitted by a single power-law spectrum with a spectral index of 
-1.59 ± 0.06. For the western relic we find a spectral index of 
-1.49 ± 0.12. We note that the 241 MHz flux measurement and 
the corresponding error are somewhat difficult to estimate as the 
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Source C. GMRT 610 MHz contour levels are displayed 
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Fig. 13. Source D. GMRT 610 MHz contour levels are displayed 
as in Fig.fTOl 



noise in the image increases sharply towards bright radio source 
E because of residual calibration errors, see also Fig. 

As reported in Sect. 13.41 the spectral indices at the front of 
the relics are about - 1 .2 ± 0.2 and - 1 .0 ± 0. 1 5 for RE and RW, re- 
spectively. The errors give the variation in spectral index across 
the outer edges of the relic. The integrated spectral indices are 
consistent with this values (if we assume that DSA takes place), 
being steeper by about 0.5 units. We take the values at the front 
of the relic as the injection spectral indices. This then gives Mach 
numbers of 2.2 + q ? and 2.4 + jj-j for relics RE and RW, respectively 
(Ivan Weeren et a i. 2009a). The uncertainties in the Mach num- 
ber are based on the variation in spectral index at the front of the 
relics of about -0.15 units. These Mach number are in agree- 
ment with those found in other merging clusters, typically be- 



ing between 1.5 and 3 (Markevit ch et alJ2 002. 20051 iMarkevitchl 
2006: lRussell et al.ll2010l:lFinoguenov et alJl2010h . 

Compression of fossil radio plasma 

(EnBlin&Gopal-Krishna 2001) does not seem to be a likely 
scenario to explain the relics in ZwCl 0008.8+5215, because 
regions with a spectral index of -0.9 are seen in front of the 
relic, not very steep for a radio phoenix. The integrated radio 
spectra do also not reveal significant spectral curvature (see 
for example fsiee et al.ll200ll) . Furthermore, relic RE has a size 
of 1.4 Mpc and the time to compress such a large radio ghost 
would have removed mo st of the energetic partic les responsible 
for the radio emission (iClarke & Ensslinl [20061) . Relic RW is 
located not far from the complex sources E and F. If RW is 
directly associated with E and F, one would not expect the 
spectral index to steepen across RW in the direction of E and 
F. Also, the overall polarization fraction of RW is considerably 
higher than that of sources E and F. Therefore, our preferred 
scenario is that of relics tracing shock fronts where particles are 
accelerated or re-accelerated by the DSA mechanism. 

The large difference of a factor of five in linear extent for 
relic RE and RW, is unlike that of previously known double 
relics. This large ratio in linear size could be explained by a rel- 
atively large mass ratio between the merging clusters, as the size 
of the shock waves formed during a binary cluster merger event 
roughly scale with the radii/masses of the merging components. 
The galaxy iso-density contours in Fig. [2] show two subclusters, 
with the eastern subcluster being slightly larger. Another possi- 
bility is that the shock front on the west side of the cluster is 
partly broken up due to the presence of substructures or galaxy 
filaments. In case of re-acceleration, a supply of fossil electrons 
is needed. Therefore, a third possibility is that if these electrons 
have a limited spatial distribution, the shock front might be illu- 
minated only locally. 

The overall configuration of the relics is largely symmetric 
around the east-west merger axis, indi cating that the imp act 



t pa - 

m. 



rameter for the merger is close to zero (IRicker & Sarazin 2001 



4.2. Radio luminosity profile for the eastern relic 

The width of a relic (/ re ii c ), tracing a plane shock wave, to first 
order reflects the downstream velocity of the shock V2, a charac- 
teristic timescale due to spectral ageing Ageing, and the angle 
between the shock front normal and the plane of the sky 



4elic 



v 2 X U 



gemg 



COS0 



(2) 



The characteristic timescale due to spectral ageing is given by 



'ageing 



[yr] = 3.2 x 10 



10 



[(l+z)v] 



-1/2 



(3) 



with B the magnetic field at the shock front in //Gauss, Bqmb 
the equivalent magnetic field strength of the CMB in //Gauss, 
and v the observed frequency in MHz. At z = 0.103, Bqmb is 
4.0 //Gauss. If V2 and cp are known this gives a method for deter- 
mining the magnetic field strength. Even if is not known, limits 
on the magnetic field can be obtained if the observed width (/ re ii c ) 
is smaller than the maximum width allowed from Eq. [2] 

To get an estimate of V2, we use a temperature in the post- 
shock region of 6 keV, i.e, about twice the average cluster 
temperature. This factor of two increase is roug hly what has 
been observed in othe r clusters with shocks (e.g.. iRussell et all 
2010t lMaet al.ll2009l). We use Rankine-Hugoniot jump condi- 
tions (lLandau & Lifshitzl 19591) . with an adiabatic index y = 5/3, 
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Fig. 15. Integrated fluxes for the relics. We find a spectral index 
of -1.59 for relic RE and -1.49 for relic RW by fitting a single 
power-law spectrum through the flux measurements at 241, 610, 
1382, and 1714 MHz. 
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Fig. 16. Radio luminosity profile across the width of relic RE. 
Solid line displays the 1382 MHz luminosity profile, averaged 
over the full 1 .4 Mpc extent, of the relic RE. Dashed line displays 
the resulting spectral index profile (labeling on the right axis), 
also averaged over the full 1.4 Mpc extent, between 1382 and 
241 MHz. The averaging was done by adding up the total flux 
at each frequency in a spherical shell and then calculating the 
spectral index. 

and take the Mach number from the injection spectral index. This 
gives 

J2 _ 5M 4 + 14M 2 - 3 

¥ x ~ 16Af ' ( } 

with indices 1 and 2 referring to the pre-shock and post-shock 
regions. The downstream speed is given by V2 = Mc St \/C, with 
c St \ the pre-shock sound speed, (7&b^iMh/^) 1/2 \ with \i - 0.6 
the mean molecular weight. The compression ratio C is given by 

13 1 

C AM 2 + 4 ' ( } 



Filling in the numbers gives C = 2.4 and c St \ - 1100 km s . 
We then obtain V2 = 750 km s" 1 . The downstream velocity 
depends only weakly on the adopted downstream temperature. 
For example, using a downstream temperature of 10 keV in- 
creases of V2 to about 950 km s" 1 . For the remainder we will 
adopt V2 = 750 km s" 1 . This then gives for the width of the relic 
(FWHM) observed at 1382 MHz 

B l l 2 , 

/relic [kpc] « 628 X — COS , (6) 

B +5 CMB 

with the magnetic field strengths in units of //Gauss. 

For = 0°, the maximum width is 46 kpc, which corre- 
sponds to B « 2 //Gauss. This is smaller than the observed width 
of about 150 kpc (see Fig. [16]) and hence no constraints on the 
magnetic field can be put since the angle is not known. It is 
possi ble to set limits on using the observed polarization frac- 
tion (lEnsslin et al.| [l998). A 20% polarization fraction implies 
< 50°. This limit on is not consistent with the observed 
width which would require > 72°. Although, for large sec- 
tions of the relic the polarization fraction is unknown and could 
be smaller than 20%. 

4.3. Simulated radio luminosity and spectral index profiles 

In the above analysis we assumed that a relic traces a planar 
shock wave. In a more realistic model of a relic would trace a 
shock wave that forms a part of a sphere. This is illustrated by 
the curved shape of relic RE. The observed width is about a fac- 
tor of three larger than the maximum intrinsic allowed width, 
max(/ re ii c (0 = 0, B)). This implies that projection effects proba- 
bly play an important role. The questions is then why do we still 
see a clear spectral index gradient (Fig. [I6j across the relic? 

To answer this questions we use a more realistic model of a 
shock front. The spherical shock subtends an angle *F into the 
plane of the sky and has a radius of curvature ^ pro j e cted- The total 
angle subtended by the relic is 2 V F. We compute the radio lumi- 
nosity profiles at the observed frequencies of 241 and 1382 MHz. 
The injection spectral index is taken to be -1.0. Synchrotron 
cooling processes, based on the distance of the emitting radio 
plasma from the front of the shock, which in turn depend on 
the downstream velocity V2 = 750 km s" 1 , are taken into ac- 
count. For the magnetic field we assume B = 2 yuGauss, which 
maximizes the intrinsic width of the relic to 46 kpc. A spectral 
index profile is computed using the profiles at the two different 
frequencies. The resulting intrinsic luminosity profiles (with in- 
trinsic referring to a planar shock wave without any projection 
effects) and profiles for 7? pro jected = 0-75 and 1.0 Mpc, with open- 
ing angles *F = 22, 30, 40° are shown in Figs. [T7] and Q21 For 
^projected = 0.75 Mpc, we find that the profile with an opening 
angle between 30 and 22°(~ 26°) provides the best match to the 
observed profile. For ^ pro j ecte( j = 1.0 Mpc, we find the best match 
for v F = 22°. 

Our computed luminosity profiles do no provide a very good 
match to the observed profile at distances of more than 0.85 Mpc 
from the cluster center. The observed profile is more symmet- 
ric, while the computed profiles are rather asymmetric with a 
strong luminosity decrease at large radii. This may be caused 
by the fact that the actual 3D shape of the shock front differs 
somewhat from a sphere. Also, we assumed a uniform surface 
brightness over the front of the shock surface (which forms part 
of a segment of a sphere). At the edges of this surface the radio 
emission drops to zero abruptly. This causes the discontinuities 
in the modeled profiles inwards of the peak luminosity towards 
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Fig. 17. Radio luminosity profile for 7? pro j ecte( j = 0.75 Mpc and 
B = 2 yuGauss. Solid black line shows the observed luminos- 
ity profile for relic RE at 1382 MHz, while the dashed black line 
represents the spectral index profile between 241 and 1382 MHz. 
The simulated intrinsic luminosity and spectral index profiles 
(i.e., without any projection effects), are shown in purple by solid 
and dashed lines, respectively. The simulated profiles, including 
the projection effects, are shown by solid lines for opening an- 
gles of 22°(cyan), 30°(green), and 40°(red). The corresponding 
spectral index profiles are shown by the dashed lines. 



the cluster center. In the GMRT and WSRT images, the relic's 
surface brightness fades towards the northern and southern ends. 
This could (partly) be explained by the spherical shell model we 
use for the relic, as the relic's extent into the plane of the sky 
decreases at the northern and southern ends. It is also possible 
that the surface brightness across the shell decreases towards the 
edges. This effect is not included in our model. However, our 
goal was not to reproduce the exact profile of the relic, but rather 
to show that although projection effects can be significant, a clear 
spectral index gradient can remain. Even for an opening angle of 
40° (a total of 80° into and out of the plane of the sky), a steep- 
ening of more than 0.5 units in the spectral index is predicted 
towards the cluster center. 

Based on this we argue that although relic RE is widened 
significantly by projection effects, the fact the we see a clear 
spectral index gradient is not surprising. This could also ex- 
plai n the spectral index gr a dients visible for the reli c s observed 
by Rottgering ir5]dl997b: IClarke & Ens^(l2006[):lQrru et al.l 
(2007h: iGiacintucci etal.1 (120081) : iBonafede et al.l (120091) . even 
though the observed widths are significantly larger than the max- 
imum intrinsic widths. 



4.4. Equipartition magnetic field strength 

Since the width of the relic is larger than the maximum intrinsic 
width, we estimate the magnetic field at the location of the relics 
by assuming minimum energy den sities in the relic s . We us e the 
same procedure as described in Ivan Weeren et al.l (l2009al) and 
take k = 100, i.e, the ratio of energy in relativistic protons to that 
in electrons. For relic RW, we have a spectral index of -1 .49, and 
a surface brightness of 1 .2 //Jy arcsec" 2 . We take 290 kpc for the 
depth (d) along the line of sight. This gives B Qq = 3.4 //Gauss. 
For RE we have a spectral index of -1.59, a surface brightness 
of 0.96 //Jy arcsec -2 , and we assume d = 1 Mpc. This gives 



Fig. 18. Radio luminosity profile for /^projected = 1.0 Mpc. See the 
caption of Fig. [T7] for an explanation of the various profiles. In 
addition, luminosity and spectral index profiles for an opening 
angle of 15° are shown in blue. 



B eq = 2.5 //Gauss. The equipartition magnetic field strength 
scales with (1 +k) 2n . In the above calculation, we used fixed fre- 
quency cutoffs (Vmin = 10 MHz and v max = 100 GHz), which i s 
not entirely correct dBeck & Krausdl2005t iBrunetti et alll 19971) . 
With low and high energy cutoffs (y min , y max ), y min <<c y max , 
and fixing y m ] n to 100, we find a revised magnetic field strength 
(B' eq ) of 7.9 and 6.6 //Gauss for RW and RE, respectively. For a 
lower cutoff of y m [ n = 5000, we get 1.4 and 1.0 //Gauss for RW 
and RE, respectively. The revised equipartition magnetic field 
strength (B' eq ) scales with (1 + &) 1/(3_a) , for different values of k. 



5. Conclusions 

We discovered a double radio relic in the galaxy cluster 
ZwCl 0008.8+5215, located at z = 0.103 (based on a single 
spectroscopic redshift). A ROSAT X-ray image and galaxy iso- 
density map show that the cluster is undergoing a binary merger 
event, with the merger axis oriented roughly east- west. The two 
radio relics are located along this merger axis, while their ori- 
entation is perpendicular to this axis. The relics probably trace 
shocks waves in the ICM, created by the merger event, in which 
particles are (re)accelerated by the DSA mechanism. Integrated 
radio spectra are consistent with particle acceleration in the 
shock by DSA and indicate Mach numbers of ~ 2 for the shocks. 
The spectral index for both relics shows a steepening towards the 
cluster center. Parts of the relics have a polarization fraction in 
the range of 5 - 25%, but further observations are needed to bet- 
ter map the polarization properties. The relics have an extent of 
1 .4 Mpc and 290 kpc. This factor of five difference in their linear 
extent is unlike that of other known double relic systems. The 
size difference could be related to a relatively large mass ratio 
between the two merging clusters, although galaxy iso-density 
contours do not indicate a large difference in masses between the 
two subclusters. Alternatively, the second shock front is partly 
broken up due to interaction with substructures, or the small size 
of the western relic reflects the limited spatial distribution of fos- 
sil electrons. 

We modeled the radio luminosity and spectral index profiles 
of the eastern relic, assuming that the relic traces a curved shock 
front. We conclude that projection effects play an important role 
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in increasing the observed s width of the relic. However, we 
find that a clear spectral index gradient remains visible for large 
opening angles. 

Future X-ray observations will be needed to investigate the 
dynamical state of the cluster, determine the mass ratio of the 
merging systems, and search for shock waves associated with 
the relics. 
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